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Abstract 

We  provide  a  new  way  of  understanding  Golay  pairs  (of 
length  N)  of  sequences  in  terms  of  the  (2N  +  l)-dimensional 
discrete  Heisenberg- Weyl  group  over  the  field  Z2.  Our 
methodology  provides  a  different  insight  into  the  nature 
of  these  sequences,  as  well  as  a  mechanism  for  designing 
sequences  with  desirable  correlation  properties.  Libraries 
of  waveforms  formed  using  these  constructions  are  able  to 
provide  collections  of  ambiguity  functions  that  cover  the 
range-Doppler  plane  in  an  efficient  way,  and  thus  provide 
the  basis  for  a  suite  of  waveforms  optimized  for  extraction 
of  information  from  the  environment  in  an  active  sensing 
context. 

1  Introduction 

Golay  complementary  sequences  are  pairs  of  sequences  of 
unimodular  complex  numbers,  usually  real,  with  the  prop¬ 
erty  that  the  sum  of  their  individual  auto-correlation  func¬ 
tions  forms  a  delta  spike  or  thumb  tack.  They  were  dis¬ 
covered  independently  by  Golay  [7],  and  by  Shapiro  [16]. 
These  sequences  and  variants  of  them  have  been  considered 
for  use  by  many  authors  in  the  construction  of  phase-coded 
radar  waveforms  and  in  modulation  schemes  for  commu¬ 
nications  (see,  for  example,  [9,  8,  15,  14,  11]).  They  have 
been  used  and  investigated  in  numerous  situations  both  in 
applications  and  in  purely  theoretical  investigations.  Effec¬ 
tively,  in  a  radar  context  they  produce  zero  range  sidelobes 
at  zero  Doppler. 

Our  aim  is  to  introduce  a  completely  new  framework  for 
understanding  Golay  pairs  (of  length  2”*)  in  terms  of  the 
(2m-|-l)-dimensional  discrete  Heisenberg- Weyl  group  over 
the  field  Z2.  This  group  is  more  familiar  in  coding  theory 
for  communications,  where  it  has  been  used  by  Calderbank 
et.  al  [4]  in  the  study  of  Kerdock  and  Preparata  codes.  The 
Heisenberg- Weyl  group  over  Z2  provides  the  mathematical 
framework  for  construction  and  analysis  of  first  and  second 
order  Reed  Muller  codes,  whereas  the  more  complicated 
Heisenberg- Weyl  group  over  Z2m  provides  the  mathemati¬ 
cal  framework  for  the  description  of  Golay  pairs.  The  latter 
group  is  the  natural  finite  version  of  the  continuous  Weyl- 
Heisenberg  group,  which  is  the  basic  object  for  the  alge¬ 


braic  treatment  of  the  radar  ambiguity  function  (see  [13]). 
We  have  been  able  to  show  how  to  pass  between  the  al¬ 
gebras  of  these  two  discrete  groups,  thus  allowing  us  to 
formulate  the  Golay  property  in  the  simpler  group  and  to 
analyse  the  structure  of  Golay  pairs  at  this  level.  This  trans¬ 
fer  property  is  encapsulated  in  an  apparently  novel  formula 
that  relates  translation  in  the  Z2  world  with  translation  in 
the  Z2™  world,  in  the  context  of  Heisenberg- Weyl  groups. 

Using  this  analysis  we  are  able  to  discuss  not  only  com¬ 
plementary  pairs,  but  also  collections  of  them  that  relate 
naturally  to  each  other.  These  collections  have  also  been 
discovered  previously  in  the  work  of  Welti  [17]  and  in  the 
so-called  Prometheus  Orthonormal  Set  (PONS)  of  Byrnes 
[3].  Variants  of  them  obtained  by  “special”  permutations  of 
the  Welti-PONS  codes  have  been  described  and  used  in  a 
communications  context  by  Popovic  [15]  and  Budisin  [2]. 
These  variants  share  the  correlation  properties  of  the  origi¬ 
nal  Welti-PONS  codes.  They  form  special  subgroups  of  the 
Z2 -Heisenberg- Weyl  group  of  operators. 

Our  work  shows  how  all  of  these  collections  arise  naturally 
from  the  algebra  of  the  situation  governed  by  the  Z2  view, 
even  though  the  Golay  property  is  describable  in  the  1^2-^ 
world.  We  have  been  able  to  characterize  the  waveform 
collections  encompassed  by  the  Popovic-Budisin  construc¬ 
tion,  applied  to  the  Welti-PONS  waveform  collections,  thus 
showing  their  special  nature  among  waveforms/codes.  Our 
methodology  provides  new  insight  into  the  nature  of  these 
codes,  and  a  new  technique  for  their  analysis,  as  well  as  a 
mechanism  for  designing  sequences  with  desirable  correla¬ 
tion  properties  both  for  communications  and  sensing  appli¬ 
cations. 

The  paper  is  organised  as  follows.  We  begin  by  de¬ 
scribing  the  general  construction  and  properties  of  finite 
Heisenberg- Weyl  groups.  In  Sections  3  and  4  we  show 
how  the  conditions  for  a  pair  of  sequences  to  be  Golay 
complementary  can  be  formulated  in  terms  of  two  partic¬ 
ular  Heisenberg- Weyl  groups.  In  Section  5  we  give  suf¬ 
ficient  conditions  for  orthonormal  bases  of  cyclic  or  Z2m- 
Golay  complementary  sequences  to  exist  and  then  demon¬ 
strate  how  to  relate  this  condition  back  to  the  usual  linear 
or  Z-Golay  complementary  sequences  in  Section  6.  As  an 
example,  in  Section  7,  we  show  how  the  Budisin-Popovic 
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Golay  complementary  sequences  [2]  fit  into  our  scheme. 

Finally,  we  note  that  this  paper  is  really  a  summary  of  re¬ 
sults.  We  have  omitted  many  of  the  proofs,  and  derivations 
have  been  shortened  or  left  out  entirely. 

2  Discrete  Heisenberg- Weyl  Groups 

In  this  section  we  give  a  short  summary  of  the  construc¬ 
tion  and  properties  of  finite  Heisenberg- Weyl  groups.  For  a 
more  detailed  description  and  further  references  see  [10]. 

We  begin  by  defining  a  configuration  space  A  =  Z™  con¬ 
sisting  of  m-tuples  of  elements  from  the  integers  modulo  g. 
In  this  paper,  we  will  take  q  =  p”,  for  some  prime  number 
p.  Under  addition,  A  forms  an  Abelian  group.  In  radar  the¬ 
ory  the  space  A,  with  m  =  1  would  represent  discretized 
version  of  the  range  space,  while  in  discrete  quantum  me¬ 
chanics  the  space  A  could  represent  possible  discrete  posi¬ 
tions  for  a  particle. 

Define  a  Hilbert  space  H,  having  orthonormal  basis 

{|a)  :  a  €  A},  (1) 

which  we  refer  to  as  the  Dirac  basis.  Note  that  we  use  the 
“bra-ket”  notation  for  elements  of  the  Hilbert  space.  An 
arbitrary  element  \(j>)  ^Ti  can  be  expanded  in  this  basis  as 

k)  =  XI  (2) 

aeA 

where  (•  |  •)  is  the  inner  product  on  Ti. 

The  dual  group  of  A,  denoted  A,  consists  of  the  homomor- 
phisms  from  the  group  A  into  the  unit  circle  H  in  C.  A  is 
also  an  Abelian  group  (under  multiplication)  and  is  (since 
A  is  finite)  isomorphic  to  A.  This  isomorphism  is  made  ex¬ 
plicit  through  identification  of  each  b  G  A  with  a  7b  G  A, 
such  that 

7b(a)  =  (3) 

for  all  a  G  A,  where  co  =  exp(27rf/p),  and  •  denotes  the 
usual  dot  product  on  Z™.  We  see  from  (3)  that  the  elements 
of  A  are  just  discrete  sinusoids,  or  multi-dimensional  ver¬ 
sions  of  such.  To  each  element  of  7b  G  A  we  can  assign  a 
vector  in  H  by 

The  set  { I  a)  :  a  G  A}  also  forms  an  orthonormal  basis  for 
Ti.,  which  we  refer  to  as  the  Fourier  basis.  We  can  define  the 
unitary  Fourier  transform  operator  relating  this  orthonormal 
basis  to  (1)  by 

jm/2 

^  o;‘’-|a)(b|,  (5) 


where  |  a)  (b  |  represents  the  cross  projection  operator  on  H 
whose  action  on  |  ^)  G  is  |  a)  (b  |  |  (/))  =  (b  |  <())  |  a) 

We  will  denote  the  group  A  x  A  ~  A  x  A,  which  we  regard 
as  a  Zq-module  over  the  ring  Z^,  by  E.  When  q  is  prime, 
Zq  is  a  field  and  is  a  vector  space.  We  will  refer  to  E  as 
the  phase  space. 

On  Fi  we  define  the  unitary  operators  {D(a,  b)  :  (a,  b)  G 
E}  by 

£>(a,b)  =  ^  |c-t-a)(c|.  (6) 

ceA 

Two  such  operators  have  the  multiplication  rule 

i:)(a,b)i:)(a',b')  =u;‘’  ‘"'i:'(aH-a',b-f  b'),  (7) 

from  which  we  have  the  commutator  relation 

D{a,h)^D{a\h')D{a,h)D{a',hy  =  I ,  (8) 

where  ^  denotes  adjoint. 

The  set  of  unitary  operators  on  Fi 

E  =  {T(A,  a,  b)  =  uj^D{a,  b)  :  A  G  Z„  (a,  b)  G  E}, 

(9) 

if  9  2,  or 

E  =  {r(A,  a,  b)  =  i^D{a,  b)  :  A  G  Z4,  (a,  b)  G  E}, 

(10) 

if  <7  =  2,  forms  an  representation  of  the  discrete 

Heisenberg- Weyl  group  on  FI.  This  representation  is  ir¬ 
reducible  [12,  4].  This  means  that  there  are  no  nontrivial 
subspaces  of  Ft  invariant  under  the  action  of  E. 

The  center  of  the  group  E,  Z{E),  consists  of  the  elements 
{to^I :  A  G  Zg},  or  {i^I :  A  G  Z4}  if  q  =  2,  where  /  is  the 
identity  operator  on  FC.  The  factor  space  E/Z{E)  is  easily 
seen  to  be  the  phase  space  E.  Considering  the  commutation 
relation  (8),  we  can  define  the  symplectic  inner  product 

((a,b),(a',b'))  =a-b'-a'-b,  (11) 

on  the  phase  space  E,  and  note  that  two  operators  D{a,  b) 
and  D{a' ,  b')  commute  if  and  only  if  ((a,  b),  (a',  b'))  =  0. 
We  may  then  identify  subgroups  of  E  consisting  mutually 
commuting  sets  of  operators  Z)(a,  b)  with  isotropic  sub- 
modules  of  E.  We  also  define  the  symplectic  dual,  or  just 
dual,  of  any  submodule  (subspace)  El  C  to  be 

=  {(a,b)  G  E  :  ((a, b),  (a',  b'))  =  0,V  (a',b')  G  H}. 

_  (12) 

An  isotropic  subspace  satisfies  H  C  H  .  It  is  maximal 
isotropic  if  and  only  if  this  inclusion  is  an  equality.  An 
isotropic  submodule  (subspace)  H  of  E  corresponds  to  the 
Abelian  subgroup  {I?(a,  b)  :  (a,  b)  G  H}  of  E. 

Finally,  in  this  section,  we  consider  the  space  of  linear  op¬ 
erators  O  on  the  Hilbert  space  Ft.  We  have  the  following 
theorem,  which  can  be  prove  by  substituting  (6)  into  (13): 


Theorem  1.  Any  operator  S  G  O  can  be  represented  as 

^=4  E  Trp(a,b)t5)i?(a,b).  (13) 

(a.b)e® 

We  refer  to  s(a, b)  =  Tr(il>(a, b)^S')  as  the  Weyl  trans¬ 
form  of  S. 

Equation  (13)  implies  that  the  map  S  — >  s(a,  b)  = 
Tr{D{a,hy S)  gives  an  isometry  from  O  to  L^(£^),  with 
the  inner  products  related  by 

Tr(5''fi?)  =  1^  s(a,b)r(a,b).  (14) 

(a,b)e£ 

We  note  that  the  orthogonal  projection  on  to  a  vector  |  (j)), 
namely  has  the  expansion 

=  ^  A4a,h)D{a,h),  (15) 

(a,b)e® 

where  ^^(a,  b)  =  {(f)\D{a,h)\(f))  is  the  ambiguity  func¬ 
tion  of  I  ^) . 

3  Golay  Sequences 

Consider  two  unimodular  sequences  of  complex  numbers 
X  and  y  of  length  N.  Two  such  sequences  are  said  to  be 
Golay  complementary  if  the  sum  of  their  respective  auto¬ 
correlation  functions  satisfy 

corrfc(x)  -f  corrfc(y)  =  2NSk,o,  (16) 

forfc  =  — (A^— 1),  •  •  •  ,  (iV— 1).  We  note  that  the  remaining 
shifts  are  automatically  zero.  Such  sequences  have  an  ex¬ 
tensive  literature,  a  sample  of  which  are  [17, 16,  8,  2, 6,  7,  1, 
5].  We  can  write  the  condition  (16),  at  least  for  cyclic  con¬ 
volutions,  in  terms  of  a  particular  finite  Heisenberg- Weyl 
group.  We  will  show  how  results  for  linear  convolutions 
can  be  recovered  in  Section  6. 

We  begin  by  considering  the  case  in  which  the  configura¬ 
tion  space,  which  label  positions  in  the  sequence,  takes  the 
form  A  =  Z2™.  In  this  case  we  define  the  Hilbert  space 
T-Cm  with  Dirac  basis 

Ad  =  {\j):J  (17) 

We  can  then  define  the  finite  Heisenberg- Weyl  group 

Sni  =  {uJ^A{j,k)  ■  M  G  ^2",  {j,k)  G  £m},  (18) 

where  Em  =  ^2"*  x  ^2™^  ^  =  exp(27ri/2’”)  and  the  uni¬ 
tary  operators  A(j,  k)  are  defined  by 

A(j,fc)=  ^  (19) 


As  noted  above  this  is  an  irreducible  representation. 

The  conditions  for  two  vectors  |  (/)),  1 7/>)  G  Tim  to  be  Z2™- 
Golay  complementary  are, 

1 .  I  (/>)  and  I  ^/>)  have  unimodular  coefficients  (up  to  an 
overall  normalising  constant),  and 

2.  for  all  non-zero  j  e  Z2m , 

(</)|A(j,0)|<^)  +  (t^|A(j,0)|V')  =  0.  (20) 

The  condition  (20)  can  be  written  equivalently  as, 

Tr((|(/.)(((.|  +  |V')(V'l)A(j,0))  =  0,  (21) 

which  implies  that  the  projector  Q  =  \  (t)){(t)\  +  \i^){4’  \  must 
be  orthogonal  to  the  subspace  5a  spanned  by  the  orthonor¬ 
mal  set  {A(j,  0)  :  j  G  ^  0}  in  the  Hilbert-Schmidt 

class  O  of  operators  on  Tim- 

Suppose,  now,  that  we  have  bijection  l  :  Z2™  ^  Z™.  Then 
each  I  j)  G  To  can  be  labelled  by  a  unique  a  G  Z2m;  that 
is,  a  binary  string,  so  that  |  j)  =  |a),  with  t(j)  =  a.  We 
take  L  to  be  the  mapping  from  j  to  its  binary  representation. 
The  bijection  l  now  allows  us  to  define  a  unitary  irreducible 
representation  of  the  finite  Heiseberg-Weyl  group  with  con¬ 
figuration  space  Z™  on  Tim'- 

Em  =  Dia,  b)  :  A  G  Z4,  (a,  b)  G  Em},  (22) 

where  Em  =  x  Z™  and  the  unitary  operators  iA(a,  b) 
are  defined  by 

77(a,b)  =  ^  (-l)‘’''^|c0a)(c|.  (23) 

Here  we  denote  addition  on  Z™  by  0  to  distinguish  it  from 
addition  on  1^2^  in  (18).  Obviously,  the  bijection  i  induces 
a  bijection  L2  ■  Em  Em  given  by  t2(a,  b)  =  {j,  k). 

We  shall  consider  the  expansion  of  the  A(j,  0)  in  terms  of 
the  79  (a,  b).  Let  us  write  this  as 

A(j,0)=  Y  aj(a,b)i9(a,b),  (24) 

(a,b)e®™ 

where  the  aj  will  be  determined  later. 

In  term  of  the  Weyl  representation  (21)  becomes 

Y  (Al0  0Alv,)(a,b)  aj(a,b)  =  0,  jV  0,  (25) 

(a,b)eE™ 

where  A(j,  and  are  the  ambiguity  functions  of  |  (/>)  and 
I  tp),  and  we  have  used  (15).  We  see  that  a  sufficient  con¬ 
dition  for  the  pair  (|  </>),  \ip))  to  be  Z2m -Golay  complemen¬ 
tary  is  that  the  supports  of  and  aj  be  disjoint 


for  all  j  ^  0.  We  will  define  the  support  of  an  opera-  4  The  Support  of  Sa 
tor  S'  G  O  to  be  the  support  of  its  Weyl  transform  s,  i.e., 

supp  S  =  supp  s,  where  Let  Sa  C  O  be  the  subspace  spanned  by  the  operators 


s(a,b)  =  ^Tr(Sf?(a,b)t),  (26) 

and  we  will  further  define  the  support  of  a  subspace  S  C  O 
to  be 

suppS  =  suppS.  (27) 

ses 

Proposition  1.  A  sufficient  condition  for  a  pair  of  vectors 
(If),  If))  to  be  ’L'f-Golay  complementary  is  that 

supp  {A4,  +  A^p)  n  suppSa  =  0.  (28) 

We  note  that  the  projector  -D(a,  b)(5_D(a,  b)^  has  the  same 
support  as  Q,  for  all  (a, b)  G  Em-  Thus,  if  {\4>),\f’)) 
satisfies  (28),  then  so  do  (-D(a,  b)|  (/)),  D{a.,  b)  1 1/;)),  for  all 
(a,  b)  G  Em- 

Given  the  last  statement  it  makes  sense  to  try  to  minimise 
the  support  of  the  projection  Q  =  |  </>)((/)  |  -f  |  |  or 

equivalently  the  function  Ap  +  Atp.  Now  the  1-D  projector 
\<l>){4>\  associated  with  maximal  isotropic  subspaces  of  Em 
or  equivalently  maximally  commutative  subgroups  of  Em, 
have  the  smallest  support  of  any  of  1-D  projectors.  In  fact, 
each  maximal  isotropic  subspace  H  C  Em  is  associated 
with  an  orthonormal  basis, 

E  =  =  D{a,h)\f)  :  (a,  b)  G  (29) 

where  |  f)  is  any  normalised  vector  in  Hm  for  which 

D{a,h)\f){f\D{a,hy  =  \f){f\,  (30) 


{A(j,0):  jGZ2™,j^0}.  (34) 

The  operator  A(j,  0)  has  expansion  coefficients  in  terms  of 
the  operators  D{a.,  b),  given  by 

a,(a,b)  =  2-™Tr(A(j,0)77(a,b)t)) 

=  2-'"  E 

kG'^2'^ 

=  2-™  E  (-1)^’'^"®®  (c  I  c  0  a  +  j) 

=  2-™  E 


We  denote  by  Ca  C  Z™  the  subspace  of  elements  of 
Z™  which  are  covered  by  a,  by  which  we  mean  that  they 
have  1  entries  only  where  a  has  1  entries.  Write  w  = 
(1,  0, 0,  •  •  •  ,  0)  G  Z^  and  w  =  (0, 1, 1,  •  •  •  ,  1)  G  Z^. 
For  j  G  Z2™ ,  let  Aj  be  the  set  of  a  G  Z™  such  that 

c  —  c  0  a  =  j,  (36) 


has  a  solution  c  G  Z™. 

Proposition  2. 


A(j,0)  = 


E 

ae^j 


1 

ICa  n  Cwl 


beCanC^i^ 


(37) 


where  Co(a,  j)  is  the  unique  solution  of 


for  all  (a,  b)  G  H.  For  all  vectors  in  |  f)  G  Ejy,  the  sup¬ 
port  of  I  (j)'){4>'  I  is  El. 

By  choosing  a  pair  of  vectors  associated  with  the  same 
maximal  isotropic  subspace  H,  that  is,  taking 

Qh  =  \  +  D{ai,hi)\(l>){(l>\D{ai,hi)\  (31) 

for  some  (ai,bi)  G  Em/H,  we  ensure  that  suppQ;^  C 
supp  Ap,  =  H.  In  fact,  the  Weyl  transformation  of  Qjp  is 

Tr(Qj^77(a,b)t)  =  (10  (32) 

so  that,  writing  K  =  {(0,  0),  (ai,  bi)}, 

supp  Qjj  =  H  f]  .  (33) 

Thus,  our  goal  is  to  find  a  maximal  isotropic  subspace  H 
and  a  subspace  K  such  that  H  O  K  n  supp  Sa  =  P- 

We  next  consider  the  form  of  supp  Sa- 


Co  -  Co  0  a  =  j,  (38) 

with  Co  G  Ca  n  Cw- 

Now  return  to  the  subspace  Sa  spanned  by  the  set 
{A(j,0)  :  j  G  Z2™,  j  0}. 

We  have  the  following  result. 

Proposition  3.  The  support  of  the  subspace  Sa  is 
supp^A  =  {(a,b)  :  a  G  Z™\{0}.b  G  Ca  n  Cw}  (39) 

Proof  First  note  that  0  f.  Aj  for  j  ^  0.  For  each  a  G 
Z™/{0},  there  is  at  least  one  non-zero  j  G  Z  such  that 
a  G  Aj.  For  any  such  j,  Q;j(a, b)  f  0,  if  and  only  if 
b  G  Ca  n  Cw.  □ 

Figure  1  shows  the  support  of  the  subspace  Sa  for  to  =  4 
and  TO  =  5.  We  note  that  the  support  has  the  form  of  a  pair 
of  Sierpinski  triangles. 


a 
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Figure  1:  Supports  of  the  subspace  Sa  for  m  =  4  and 
m  =  5. 


as  defined  in  (40).  If  there  exists  a  character  7(a(],bo)  '^/ 
Ziff  X  Zff,  such  that 


7(ao,bo)({(a,  Pa)  :  a  G  Z-\{0},Pa  G  Ca  n  Cw}) 

=  {-!}, 


(43) 


then  for  each  \(j))  G  IFjjp,  the  pair  (|  (/)),  D(ao,  bo)|  </>))  is 
Z^^-Golay  complementary. 

Proof.  This  follows  from  the  arguments  in  Section  3.  □ 

Note  that  if  there  exist  several  characters  7(a.b)  of  Z™  x  Z™, 
with  (a,  b)  G  Em/Hp  such  that  (43)  is  satisfied,  then  each 
I  f)  G  will  have  several  Goly  complementary  partners, 
one  corresponding  to  each  character.  Also,  the  condition 
(43)  is  equivalent  to  stating  that  we  can  find  an  operator 
D(ao,bo)  which  anticommutes  with  every  element  of  the 
intersection  {D(a,  b)  :  (a,  b)  G  Hp  n  supp^A}. 

Example  For  m  =  4  Figure  2  shows  the  supports  of  the 
maximal  isotropic  subspaces  Hp  and  Hq  with 


1 

0 

1 

1 

1 

0 

0 

0 

and  Q  = 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

Vi 

0 

0 

0) 

0 

1 

0/ 

(44) 

respectively.  The  maximal  isotropic  subspaces  are  show  in 
red  and  magenta,  with  magenta  at  the  points  of  overlap  with 
Sa-  Hq  has  Z2™-Golay  complementary  sequences  while 
Hp  does  not.  For  Hq  two  characters  satisfying  (43)  are 

7((o.o.o,o),(o.o.o.i))  7((o.o.o,o),(o,i,o.o))- 


5  Orthonormal  bases  of  Z2™-Golay  pairs 

Let  us  now  consider  the  maximal  isotropic  subspaces  Hp  G 
Em  which  have  an  associated  orthonormal  basis  consist¬ 
ing  of  elements  whose  coefficients  have  constant  modulus. 
Such  maximal  isotropic  subspaces  take  the  form 

i7p  =  {(a,Pa)  :aGZ^},  (40) 

where  P  is  any  binary  symmetric  matrix.  In  this  case  we 
have 

iFpnsupp^A  =  {(a,  Pa)  :  a  G  Z™\{0},  Pa  G  CanCw}- 

(41) 

The  characters  of  Z™  x  Z™  take  the  form 

7(a',bo(a,b)  =  (-l)‘^''‘’®-‘’',  (42) 

where  (a',b')  G  Z™  x  Z™. 

Proposition  4.  Let  Hp  G  Em  be  the  maximal  isotropic 
subspace  corresponding  to  the  binary  symmetric  matrix  P 


6  Z-Golay  Sequences 

We  are  particularly  interested  in  conventional  or  Z-Golay 
sequences,  where  the  autocorrelation  of  the  sequences  is 
calculated  on  Z  rather  than  on  Z2™ .  We  now  give  a  method 
for  constructing  such  sequences  from  sequences  which  are 
Golay  on  Z2m+i .  We  decompose  Z™®^  into  two  subspaces 
=  L  (B  U  corresponding  to  a  G  Z^®^  whose  most 
significant  digit  is  a  0  or  a  1  respectively.  Both  L  and  U 
are  isomorphic  to  Z™.  Consider  a  vector  |  fj)  G  Hm  of  the 
form 

\i’)  =  '^PeL\a).  (45) 

a^L 

If  a  pair  of  vectors  f’’))  of  this  form  are  Z2™ -Golay 

complementary  then  the  sequences  ({/3a  :  a  G  L},{/3(,  : 
a  G  L})  will  be  Golay  complementary  on  Z. 

Consider  the  ambiguity  function  of  a  vector  of  the  form 
(45).  We  have 

A-^{a,  b)  =  0,  if  a  G  U,  (46) 


then  for  each  \(j))  S  ^Hp>  the  pair  (\4i),D{a.t),hf)\(j)))  is 
Jj-Golay  complementary. 
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1000 


1111 
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We  will  refer  to  a  maximal  isotropic  subspace  for  which 
such  a  character  exists  as  Z-Golay.  We  note  here  that  [6] 
points  out  the  connection  between  Golay  sequences  and 
certain  cosets  of  the  first  order  Reed-Muller  codes  in  the 
second  order  Reed-Muller  codes.  These  cosets  are  inti¬ 
mately  connected  with  the  above  maximal  isotropic  sub¬ 
spaces  [4,  10]. 

This  situation  is  illustrated  in  Figure  3  for  the  maximal 
isotropic  subspace  Hp  <Z  Em  with 


a 


Figure  2:  For  m  =  4,  the  maximal  isotropic  subspaces  Hp 
(upper)  and  Hq  (lower)  are  shown  in  red  and  magenta,  with 
magenta  at  the  points  of  overlap  with  5a.  The  matrices  P 
and  Q  are  given  in  (44).  Hq  has  Z2m-Golay  complemen¬ 
tary  sequences  while  H p  does  not. 


and 

A^{a,h®w)  =  A^{a,h).  (47) 

Identifying  Em  with  Lx  L,we  only  need  to  ensure  that  the 
supports  of  Ac/,  +  A,p  do  not  intersect  in  L  x  L.  That  is,  we 
replace  supp  Sa  is  the  above  arguments  with  the  set 


/O  1  0  0\ 
10  10 
0  10  1 
\0  0  1  0/ 


(50) 


Two  characters  satisfying  (49)  are  7((o,o,o.o),(o,o,o,i))  ^nd 
7((o. 0,0,0), (1,0, 0,0))- 


a 


Figure  3:  The  set  TZa  (black  and  magenta)  and  the  support 
of  the  maximal  isotropic  subspace  Hp  €  E4  with  P  given 
by  (50)  (red  and  magenta).  The  intersection  between  the 
two  sets  is  shown  in  magenta 


7^A  =  {(a, b)  :  a  G  Z^\{0},  b  G  Ca},  (48) 
a  picture  of  which  is  shown  in  Figure  3  for  m  =  4. 

We  can  state  the  following  sufficient  condition  for  the  ex¬ 
istence  of  orthonormal  bases  of  Z-Golay  complementary 
sequences. 

Proposition  5.  Let  Hp  G  Em  be  the  maximal  isotropic 
subspace  corresponding  to  the  binary  symmetric  matrix  P 
as  defined  in  (40).  If  there  exists  a  character  7(ao,bo)  of 
Z™  X  Z™,  such  that 

7(ao,b„)({(a,  Pa)  :  a  G  Z™\{0},  Pa  G  Ca})  =  {-1}, 

(49) 


7  An  Example  —  Budisin  Golay  Sequences 

In  this  section  we  consider  the  maximal  isotropic  subspaces 
associated  with  binary  symmetric  matrices  of  the  form  (50), 
in  general,  matrices  with  elements  of  the  form 

[Pm] ^2,1  +  1  4“  ;  f  j  1 ‘  '  7  (51) 

wherever  the  indices  make  sense.  These  generate  the  clas¬ 
sical  Welti-PONS  Golay  complementary  sequences.  We 
write 

Tm  =  {(a,  P ^a)  :  a  G  Z™\{0},  P ^a  G  Ca}.  (52) 


For  m  =  2  we  have 


P2  = 


0  r 

1  0/  ’ 


and 


X2  —  {(a,  P 2a)  :  a  G  Z2\{0},  P2a  G  C^} 


=  {((1,1), (1,1))}- 

For  TO  =  3  we  have 


(53) 


(54) 


^0  1  0^ 
Pa  =  I  1  0  1 
,0  1  0, 


(55) 


and 


X3  =  {((1, 0, 1),  (0, 0, 0)),  ((1, 1, 1),  (1, 0, 1))}.  (56) 

For  TO  =  4  we  have 


P4  = 


/O  1  0  0\ 
10  10 
0  10  1 
\0  0  1  0/ 


(57) 


and 


X4  ={((1, 1, 0, 1),  (1, 1,0,0)),  ((1,0, 1, 1),  (0,0, 1, 1)), 
((1, 1,1,1), (1,0,0,!))}. 

(58) 


The  Budisin  Golay  complementary  sequences  correspond 
to  the  maximal  isotropic  subspaces  with  binary  symmetric 
matrices 

P'^  =  SP^S^,  (63) 

with  Pm  given  by  (50)  with  S  any  permutation  matrix. 
There  are  to!/2  distinct  P'.  If  we  add  an  arbitrary  non-zero 
binary  diagonal  matrix  to  P'm  we  also  obtain  Z-Golay  com¬ 
plementary  sequences  but  taking  values  in  {l,i,— 1,— z} 
rather  than  in  {1,  —1}. 

8  Conclusion 

We  have  show  in  this  paper  that  the  finite  Heisenberg- Weyl 
groups  provide  a  new  and  fruitful  way  of  viewing  the  Go- 
lay  complementary  sequences  (complementary  waveforms) 
and  of  understanding  the  origin  of  their  special  proper¬ 
ties.  The  work  we  have  presented  here  is  a  summary  of 
particular  results  in  a  more  general  program  of  work  in 
which  the  finite  Heisenberg- Weyl  groups  are  used  as  a  uni¬ 
fying  mathematical  structure  for  analysing  the  correlation 
and  cross-correlation  of  unimodular  sequences  and  conse¬ 
quently  phase  code  radar  pulses. 

Libraries  of  waveforms  formed  using  these  methods  are 
able  to  provide  collections  of  ambiguity  functions  that 
cover  the  range-Doppler  plane  in  an  efficient  way,  and  thus 
provide  the  basic  idea  for  a  suite  of  waveforms  optimized 
for  extraction  of  information  from  the  environment  in  an 
active  sensing  context. 


Now  we  observe  generally  that  if  the  first  or  last  element 
of  the  vector  a  is  0,  then  Pa  €  Ca  if,  and  only  if,  a  =  0. 
Furthermore,  if  a  has  a  0  at  the  1  <  rz*  <  to  element  and 
Pa  G  Ca,  then  the  (n  —  1)*  and  (n  +  1)*  elements  must 
be  1,  i.e.,  a  0  always  occurs  in  the  combination  1, 0, 1.  This 
implies  that  the  ”a”  vectors  in  Tm  are  obtained  by  append¬ 
ing  0, 1  to  the  vectors  from  Tm-2  and  appending  1  to  the 
vectors  fromZ^-i.  Thus, 

|Xm|  =  |Xm-l|  +  |Xm-2|,  (59) 

for  TO  >  3,  and  so  we  see  that  for  the  matrices  Pm  in  (51) 

|{(s,  P m^)  ■  S  G  ^™\{0},  P m^  G  C'all  =  Fm,  (60) 

where  Fm  is  the  to*  term  of  the  Fibonacci  sequence  {Fi  = 
F2  =  1).  It  is  also  easily  seen  that  for  all  to  >  1 

7((o,o.-"  ,o,o),(o,o.-"  ,o.i))(7^m)  =  {“!},  (61) 

and 

7((o,o.-"  .o,o),(i,o,-"  ,o,o))(7^m)  =  {~1}-  (62) 

So  for  all  TO  >  1  the  maximal  isotropic  subspace  Hp^  G 
Em  is  Z-Golay. 
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